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Two-scale extensions for non-periodic coefficients * 

Vsevolod Lapteva- 


Abstract 

We consider non-homogeneous media with properties which can be characterized by rapidly 
oscillated coefficients. For such coefficients we define a notion of two-scale extension, present 
several ways to construct two-scale extensions, discuss their properties and relation to homoge¬ 
nization 
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1 Introduction 

It is usually difficult to predict a global behaviour of some process in heterogeneous media (for 
example composite/porous materials) although the physics of the process might be well understood 
locally. The reason lying in the complexity of the microstructure gives rise to different upscaling 
methods. 

Heterogeneities having periodic microstructure play a central role in the development of upscaled 
models. From one side they represent an important particular case of general heterogeneous media 
and on the other there are well developed mathematical techniques (e.g. the two-scale asymptotic 
expansion method), which help to derive formally and often rigorously the upscaled model. As 
a result many physical processes in heterogeneous media having periodic microstructures are well 
investigated both from theoretical and from practical points of view and the periodicity assumption 
is usually a starting point for the upscaling procedures I2MX3MX5I- Although this assumption is 
valid in only limited number of cases, mostly in artificially created materials. Therefore for practical 
purposes one should be able to deal with non- periodic structures. 

The deterministic homogenization procedure starts from a sequence of problems {V £ }. In the 
periodic case the heterogeneity in V s is usually described by an e-periodic function a £ (x ) = a(x/e), 
where a(y) is a given V-periodic function in ( Y = (0, l) d is a period: a(y + eQ = a(y), e* is a 
unit vector, i = 1,d). Quite often the purely periodic coefficient can be generalized without 
difficulties to the locally periodic coefficient a £ (x ) = a(x,x/e) (where a(x,y) is a given V-periodic 
function in y). In the following steps one has to investigate the convergence of the sequence (in a 
wide sense) and to find a limit problem V° . The solution of the limit problem can be used in order 
to approximate the solutions of the problems V s for small enough e. 

The coefficients a(y) or a(x, y ) are considered in mathematical literature as given functions 
belonging to some functional spaces, without paying much attention where they come from. The 

‘This work was carried out during the tenures of a fellowship from University/ITWM in Kaiserslautern (Germany) 
and an ERCIM fellowship in Luxembourg and Norway. 

} NTNU-IMF, Alfred Getz vei 1, NO-7491 Trondheim, Norway (laptevv@mail.ru, vsevolod.laptev@math.ntnu.no) 


1 



construction of these coefficients which is important for usage of homogenization results will be 
discussed in this article. 

Let us assume that some process in a heterogeneous medium occupying a bounded domain 
can be described by some PDE(s) with (at least one) rapidly oscillated coefficient om{x), 
which is non necessarily periodic. This is our initial problem V. 

Asymptotical approach applied to V means that we are not going to solve it directly, but to 
construct a sequence of imaginary problems {V s } passing through V at some e: 

V £ ° ... 'P £n ~ i V £ V £n + l ... V £ ... ■> V°. 

II ’ ’ ’ ' (1) 

V 


If the sequence {V £ } is convergent in some sense to a limit problem V° which is easier than V 
then the solution of can be used to approximate (in some sense) the solutions of V £ , and in 
particular, of V £ (it is our main goal). The ’’convergence of problems” is related to convergence of 
their solutions, but it might be restrictive to say something more precise. 

In the periodic case, namely when clm(x) is e-periodic in fl there is a y-periodic function a(y) 
defined in such that clm(x ) = a(x/e). The standard sequence {V £ } is based on the e-periodic 
coefficient a e (x) = a(x/e). a £ (x) = clm{x) and consequently the condition V £ = V is not difficult 
to satisfy. This approach cannot be used for non-periodic clm(x) since there is no such periodic 
a(y) exists (except the case when the period contains the whole fl). But using the sequence {' P £ } 
based on locally periodic function a(x,y ), where the coefficients have the form a £ (x ) = a(x,x/e), 
the requirement V £ = V becomes much more realizable. We only need to find such function a(x, y ) 
and e that a(x,x/e) = clm{x). Therefore it is reasonable to make the following definition. 

Definition 1.1. Let us say that a function a(x,y), (x,y) G Ll x M d , Y-periodic in the variable y is 
a two-scale extension for om{x) if there exists a positive number e such that 


a (x, = o,m{x), Vx 


G 0. 


( 2 ) 


The article is organized as follows. In the next section several ways to construct two-scale 
extension for arbitrary initial coefficients um(x) are presented. The Section 01 contains a short 
introduction to the two-scale convergence method together with a definition and a criterion for the 
concept of admissible test function. The criterion is needed to show that the proposed in Section Q 
two scale extensions are admissible test functions in the sense of the two-scale convergence. This is 
the main purpose of Sections 0 El El (its justification consists of several results which may also be 
useful of their own). The application to the second order elliptic equation is discussed in Section [3 

Why do we need this? There are both theoretical and practical reasons to consider two-scale 
extensions. First of all, they seem to be naturally related to the formal method of two-scale 
asymptotic expansions and to its rigorous version - the two-scale convergence method. If some 
mathematical model of a physical process allows the formal homogenization procedure via two- 
scale asymptotic expansions in the case of smooth locally periodic coefficients then as the next step 
one can substitute two-scale extensions for these coefficients and check whether the homogenization 
procedure remains working for non-periodic coefficients. 

Let us now assume that our mathematical model is based on the second order elliptic equation. 
The two-scale extensions might be useful for better understanding of the following important ques¬ 
tions related to the concept of the averaged coefficient: 

- its definition, existence, properties, limits of applicability, averaging size; 
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- connection between deterministic and stochastic approaches; 

- reiterative averaging (averaging of the averaged coefficient). 

There are many algorithms currently known for practical calculation of the averaged coefficient 
(see e.g. 0,0,EH ). Some of them (having the same local problem with periodic boundary condi¬ 
tions) can be recovered by a special choice of the two-scale extension. This gives them a justification 
by an asymptotical argument as well as some freedom for improvement and generalization. For 
example it is possible to correct the averaged solution in a postprocessing step using a standard 
technique from homogenization theory [2| p.76]. Therefore for the practical problems like heat 
transfer in composite materials and unsaturated flow in heterogeneous porous media the choice of 
the two-scale extension defines a numerical method which can be used as a possible alternative to 
such methods as multiscale finite element method 0,0 or heterogeneous multiscale method 0. 

2 Three approaches to construct a two—scale extension 

First of all we have the Trivial Extension: 

a(x,y) := om(®), R d . 

But we cannot expect something better than the constant sequence {V e } = {'V } with the limit 
problem V° = V which is just as difficult to solve. This practically useless extension gives although 
an approximation to V with a perfect quality. Different two-scale extensions lead to upscaled 
problems with different quality. At least we know that not all are bad. 

For the other two approaches we need to know oj^(i) in a neighbourhood of a point in P. Since 
this can create some problems close to the boundary, let us assume that can be somehow 

extended to a larger domain P which is also bounded (if we find nothing better, we can choose 
some value of om(') hr P as a constant value in P \ P). 

Next we need to choose e. For periodic aM(x ) it is reasonable to choose e equal to the period, 
but in general we are free in choosing it. Let W{x) be an e-cube with the center x and sides 
aligned with the coordinate axes. Up to now the only restrictions on e are: we consider e to be 
small comparing to the typical size of P and all cubes W (x), i£P should be completely inside P. 

Having in mind the volume averaging method it might be reasonable to call W(x) as a (cubic) 
representative elementary volume (REV) around the point x. 

Two approaches to construct the two-scale extension a(x,y) for om(x) are different in the 
sense that the first is created via continuous ( Continuous Extension) and the second via discrete 
(" Discrete Extension) ’motion’ of W(x) in P. 

2.1 Continuous Extension 
Let x be some fixed point in P. 

• First we define an auxiliary function a(x, •) at y £ W{x): 

a(x,y) = a M (y), yeW(x). 

• Secondly we extend it to the whole R d periodically - a(x, y) is e-periodic in y. 

• Thirdly 

a(x,y) := a(x,ey) 

is defined in P x R d , V-periodic in y. It satisfies 0 . 
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2.2 Discrete Extension 

Let us assume that we have some finite partition Cl = U jClj, j = 1, ..., N. Cli n Clj = 0, i j. For 
each Qj there is a corresponding e-cube Wj = W(x 3 ), Clj C Wj. x 3 is a center of Wj. 

• First for any fixed x E Q fl Qj we define an auxiliary function a(x, •) at y E Wj-: 

a(x,y) := a M (y), y^Wj. 

• Secondly we extend it to the whole periodically - a(x, y) is e-periodic in y. 

• Thirdly 

a(z,y) := a(x,ey) 

is defined in Cl x R d , T-periodic in y. It satisfies ©. 

Remark 2.1. Both extensions are also well defined in Cl x V (this will help to show continuity of 
some properties in Cl). For Clj much smaller than Wj, x 3 E Clj the Discrete Extension can he seen 
as a discretization of the Continuous Extension. 

In order to use the results of convergence and error estimations, one usually needs smoothness 
of a(x,y). However it is easy to see that a(x,y) are continuous neither in x nor in y (and are 
properly defined only a.e.). Anyway, in the next sections our goal will be to show that these a(x, y) 
can be considered as admissible test functions in the sense of two-scale convergence and at least 
for the second order elliptic equation with highly oscillated (conductivity, permeability) coefficient 
the standard procedure T] still works and solutions of {V s } converge to the solution of V°. 

Please note that the convergence of the solutions of {V £ } is important, but it cannot guarantee 
that the solution of V can be well-approximated with the help of the solution of the problem V°. 
The approximation may fail since V £ plays a central role in the construction of the sequence and 
even if the sequence ’’converges”, V° may be ’close’ to practically useless problems V £ , for eCe 
but still ’far’ from V £ . 

2.3 An example of V s for the elliptic problem V 

In this example we consider the second order elliptic problem with homogeneous Dirichlet boundary 
condition as the initial problem 

V: - V • (a M (x)Vu) = f in fl, u\ d n = 0; (3) 

and the sequence of problems {V £ } is 

V £ : — V • (a(x, x/e)Vu £ ) = / in fl, u £ \ dn = 0, (4) 

where om( x) = {a l ^{x)\ and a(x,y) = {a l3 (x,y)} are dx d matrix functions in general case and 
a i3 (x,y ) is a two-scale extension of af( { (x). 

Naturally om(-) is required to be bounded and positive definite. Can we expect similar proper¬ 
ties for a(x, x/e) which are important for verification that {V £ } is a sequence of solvable problems? 

2.4 Properties of Continuous and Discrete Extensions inherited from om{x) 
Proposition 2.1. A property of om(x) which is valid for all x E Cl is also valid for a(x,y ) in 

n x y. 

Proof. For both Continuous Extension and Discrete Extension there is a mapping z : Cl x Y —> 0. 
that a(x,y) = aM{z(x,y)). □ 
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a(-,-). Then 


Corollary 2.1. Let Ai be the mapping Ai : om( - ) 

• At is linear. 

• |om(-)I Kv)l- 

• aikf(-) p «(•; -) p - 

• if 0-M(x) is uniformly bounded, positive definite matrix function in a^-(- ) a* J (-,-) f/ien 

a(x,y) is uniformly bounded, positive definite matrix function in Q xY. 

Proof. For example, if &m( - ) = |cim(')I then 

b(x,y ) = bM(z(x,y)) = \aM(z(x,y))\ = |a(x,y)|. Similar with others. □ 

We note that JA for the Discrete Extension has some similarity with the unfolding operator T 

0 


3 Two—scale convergence and admissible test functions 


The concept of two-scale convergence was introduced in HU and further developed in [I]. A 
recent review of a two-scale convergence in L p (Ll) space can be found in TOb In this section 
we formulate some results related to two-scale convergence in L 2 (Q) mainly following jy, but 
with some modifications of the concept of admissible test function. We will need these results in 
Section El 

Definition 3.1. Let Btf — D(fil xY) be a base space of test functions. 

A function f(x,y ) initially defined a.e. in Ll x Y we can extend to a T-periodic function in 
Ll x R. d by periodical repetition, except perhaps the points periodic to dY. 

Lemma 3.1. For any Y-periodic function ip{x,y) E C(Ll x Y) 


lim / ib(x, '—)dx 


if(x, y ) dx dy 


n jy 


(5) 


Proof. For example see 1121 . □ 

In the following we will deal with sequences { u £ }. u £ is a pair (u, e) E L 2 (fil) x M + . The sequence 
{u £ } is a sequence of pairs {(«„, e n )}^T 0 where {e n } is a fixed sequence of strictly positive numbers 
tending to zero. ”lim” is the same as ” lim ”. 

e—>0 n—KX) 

e=£n 

Definition 3.2. A sequence |u £ (x)) from L 2 (Ll) is said to be two-scale converqent to a limit 
u 0 (x,y)£L 2 (nxY) if 


(i) for all E Btf ■' 


lim / u e (x)ip(x, —)dx= / / uq(x, y)ip(x, y) dx dy (6) 

Jn £ Jn Jy 


(ii) u £ is bounded in L 2 (fil). 

We prefer to insure that all two-scale convergent sequences are bounded. Having chosen Btf 
somewhat larger, for instance T 2 [H; C per (Y)\ we would have (i) =>■ (ii) due to weak convergence of 
u £ . We refer to m for the discussion of this topic and for the definitions of the functional spaces 
like L 2 [Pt-C per {Y)]. 
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Remark 3.1. The Def. Id. 21 has sense since the limit uo(x,y) is unique as an element of L 2 (Vl x Y) 
due to density of Btf in L 2 (Vl x Y) and at least the following sequences are two-scale convergent: 

1. If cf(x,y) £ C(fl x Y) then u £ (x) = <f(x,x/e) two-scale converges to cf(x,y). 

2. If u £ (x) —* u(x) in L 2 (Vl) then u £ (x) two-scale converges to uo(x,y) = u(x). 

Proof. In both cases u £ is bounded. The first statement is a consequence of Lem. Id.II To use 
Lem. ldTl in the second statement we should approximate u(x) by a smooth function in If 2 (VI). □ 

Remark 3.2. Usually in the definition of the two-scale convergence one uses Btf = Cff r (Y)\. 

If we want to check that some sequence {u £ } is two-scale convergent then it is better to have 
possibly smaller set of test functions (Btf)- But if we already know that {u £ } is two-scale conver¬ 
gent (for example from compactness result, see Cor. 13.11) then it is desirable to be much more free 
in choosing if for ©• 

Definition 3.3. A Y-periodic function <f(x, y) square integrable inVlxY with well defined <f(x,x/e) 
in L 2 (Vl) for all e £ {£ n } is called an admissible test function (ATF) if for all two-scale convergent 
sequences {tx e } with a limit uo(x,y) holds: 


lim / u e (x)(j)(x, —) dx = / / uo(x,y)<f>(x,y) dx dy. 

Jn £ Jo. Jy 

Please note that we do not consider cf to be an element of I? (VI x Y) since different represen¬ 
tatives 4>(x,y), (f(x,y) of the same element <f £ L 2 (Vl x Y) may have 4>(x, |) 4>(x, |) in L 2 (Vl ), 

or even <j)(x, |) ^ L 2 (Vl). 

Theorem 3.1. (see m Th.1.2.) From any bounded sequence {u £ } in L 2 (VI) it is possible to extract 
a siLbsequence {?4} and there exists uo(x,y) £ L 2 (Vl x Y) so that for all 4>(x,y) £ L 2 [Vl; C per (Y)\: 


lim / u' £ (x)4>(x, —) dx = / / uo(x, y)(f>(x, y) dx dy. 

e Jn Jy 

Corollary 3.1. From any bounded sequence {u £ } in L 2 (Vl) it is possible to extract a two-scale 
convergent subsequence. 

Proof. Btf = D(Vl x Y) c L 2 [Vt-, C per (Y)]. □ 

Corollary 3.2. All functions from L 2 [Vl;Cp er (Y)] are ATF. 

Proof. Let us assume the opposite: cj>(x,y) £ L 2 [Vl-, C per (Y)\ and u £ (x) two-scale converges to 
u o(x,y), but there exists 5 > 0, subsequence {u' e } that 


/ u' £ (x)(f>(x, x/e) dx — / / uo(x,y)(f>(x,y) dx dy\ >5. 

Jn Jn Jy 


(7) 


From Def. Dii), Th.EU there exists a subsequence u £ in vJ e that for all if(x,y) £ L 2 [Vl-, C per (Y)\ 


lim / u £ (x)if(x, —) dx = / / u\(x,y)if(x,y) dxdy. 

Jn e Jn Jy 

u\(x,y) = uo(x,y) due to the uniqueness of the two-scale limit of {uf}. For if = (f> there is a 
contradiction with 0. □ 
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3.1 Necessary and sufficient conditions for 0 to be ATF 

Let us assume that (p(x, y) G Atf ~ a set of ATF. 

• First we test Def. run with u e (x) = il>(x,x/s), for all 0 G Btf■ {u £ } two scale converges to 
y) (see Rem. till) 


ip(x,y)<j>(x,y)dx dy. 


n JY 


lim [ ip(x,—)(j)(x,—)dx 
£ ^° Jn £ £ 

Second we test Def. 13.31 with u e (x) = u(x), for all u{x) G L 2 (fT) (see Rem. mi) . 

lim / u(x)cf>(x, —) dx = / / u(x)cf)(x,y) dx dy. 

J o £■ Jn Jy 

then 0(x, x/e) weakly converges to j y 0(x, y) dy in L 2 (Ll) and consequently is bounded. From 
1 and 2 we conclude that u e (x) = <fi(x,x/e) two-scale converges to cj>(x,y). 


Third we test Def. 1301 with u e {x) = <f>(x,x/e)\ 


lim / 0(x, -) 

£ ^° J n £ 


' dx = 


4>(x,y) 2 dxdy. 


n JY 


Proposition 3.1. The necessary conditions for a function <f to he from Atf : 

(j)(x,x/e) two-scale converges to <f(x,y) 


(8a) 


lim || 0 (x,x/e)|| i 2 (n) = U(x,y)\\ L 2 {nxY ) ( 8 b) 

The conditions implicitly require that 0(x, y) is square integrable in D x Y and (f>(x,x/e) G T 2 (D) 
is well-dehned for all e G {e n }- 

Theorem 3.2. Let u e {x),v £ {x) G L 2 (Q) two-scale converge to uq(x, y), vq(x, y) G L 2 (Q x Y) re¬ 
spectively. And also lim \\u £ ||z, 2 (r 2 ) = ll' u o||L 2 (nxr) then 

lim / u £ (x)v £ (x) dx = / / u 0 (x, y)v 0 (x, y) dx dy. 

£ ^°Jn Jn Jy 

Proof. See the proof of Th.1.8 in |J. There one can choose ip n (x,y) from Btf = T>(£1 x Y), 
4>{x) = 1 even if 1 is not in D(f2). Note that v £ must be bounded in assumptions of Th.1.8 and 
here it is due to (ii) in Def. l3~2l □ 

Corollary 3.3. Necessary conditions (0) are also sufficient for a function to be ATF. 

Proof. In Th. 13.21 u e (x) = 4>(x,x/e). 0 satisfies conditions ©. v £ is an arbitrary two-scale 

convergent sequence. By Def. I3.3| d> is ATF. □ 

With the help of (jHJ) we can verify whether a particular function is ATF. The condition m 
alone is not enough 1121 Rem. 1.4.5]. Although having a linear space of functions satisfying 
there is no need to check (18aK : 

Proposition 3.2. Let L be a linear space of functions such that L D Btf and all functions from 
L satisfy iSh\) . Then L c Atf- 
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Proof. We have to check (EH) for (f> E L. u e (x) = (j>(x,x/e) is bounded due to (|8b|). For any 
E B T f- 

4>(x, x/e)ip(x, x/e) = - | [</>(x, x/e) + ip(x, x/e)] 2 — </>(x, x/e) 2 — if(x, x/e) 2 j 


lim / <b(x, — )-0(x, —) dx = 

£—>o J v e' v e 

Q f2 Y 


<t>{x,y)'4’{x,y) dx dy 


We used (I8bll for f> + if, <f>, E L. 


□ 


Proposition 3.3. 4 /t_f is a linear space. 

Proof. Let <f>\. (f >2 € Atf, real numbers a, (3. We need to check © for a<f \ + /3</> 2 . For any -0 E Btf, 
EH is valid: 


lim 

£—>0 


J [{acfii + (3(f>2)ip\ (x, x/e) dx 

Q 


[(a<j>i + /30 2 )^] (®, y) dx dy. 


n y 


We can use (I8bft for <f>i , f >2 and Th. 13.21 with u e (x) = 4>i(x,x/e), v e (x) = f 2 (x,x/e) to verify EH 
for a(j)\ + /3(f> 2 - 


[a<j>i + /3</> 2 ] 2 (x,x/e) = [a 2 ^ 2 + /3 2 ^ 2 + 2a/?0i</> 2 ](x, x/e). 


Hence 


lim 

£—>0 


J [acj> 1 + /%] z 


(x,x/e)dx = / / [at^i +/?</> 2 ] 2 (x, y) dx dy. 


o y 


□ 


The following sections contain properties of the two-scale Continuous and Viscrete extensions 
of ajv/ respectively. Our main goal is to show that these extensions a(x,y) are ATF. If qm £ L 1 (H) 
then we assume that a(x,y) and a(x,x/e) are constructed pointwise a.e. in O x Y and in fI from 
some representative om{x) of a m■ Another representative <im{x) results in a.e. the same functions 
a(x,y) and a(x,x/e). 


4 Properties of the two-scale Continuous Extension 

In this section we deal only with the extension a(x,y) constructed from om(x) in the subsection 

EH1 

Proposition 4.1. For fixed x E P, a(x, ■) was constructed piecewise from aM(-), namely M rf is 
divided into \ d -cubes, by the grid 

M x {x) = jy E R d | 3k E {1,... , d}, i E Z : y k = x k /e + i - 1/2^ , 
each cube corresponds to the same e-cube W(x). 

Proposition 4.2. Let us now fix some y E R d . The function a(-,y) is piecewise constant on x E LI: 
for each y, O is divided by cubic e grid 

My{y) = {x E LI I y E A f x {x)} = {x E LI \ 3k, i : x k = y k £ - {i - 1/2)^} 


into parts where a(-,y) is constant. 


The way in which a(x,y ) was constructed makes it difficult to deal with a(x, x/e). We need 
a simple representation of a(x,x/e) for £ / £. The first argument x determines the set N x (x) 
in The second argument x/e determines which value of om(x) in the neighbourhood W(x) 
should be taken as the value a(x,x/e). The non-periodicity of om(x’) causes an uncertainty when 
x/e G A f x (x). 

Af = jx G M. d | x/e G A/Z/x) j = |x € | z : x^/e = x^/e + z — 1/2 j 

Af = jx G | 3k G {1, ... , d},i G Z : Xfc = (z — 1/2 )ee/ (e — e) j 
Af divides into open cubes A/ with a side A = ee/[e — e| and centers in 

xi = -fZf 1 ' I = (h... i d ) e Z d 

Let J e be a set of multiindexes / G that A j := Aj fl is not an empty set. 

If x/ G fl then 


a fiy, = a jxy, ^x 


a(x/,x/ + e/) = {e-periodicity} = a(x/,x/) = glm(%i ) 


Similar if x G A/ then x = x/ + h, \hk\ < A/2 


a ( x/ + h, ———j = a(xj + h, ~{xj + /i)) = {e-periodicity} = a(xy + h, xy -f—/z) = 


= clm(xi 4— /z) since x/ 4 — h G W{xi + ft,). 

We have proved the following 

Proposition 4.3. T/ze simple representation of a(x, x/e) /or all e > 0, e/e is: 

X, —J = Om{xi -|—(x — X/)) 

or using the Heaviside function 1a, (x) freing 1 in A/ and 0 elsewhere we have: 

a(x, x/e) = ^ 1 A/ (x)ai\i{xi + -(x - x/)) (9) 

/eJe e 


Roughly speaking for e < e [e > e] a(x,x/e) is built from compressed [stretched] cubes taken 
from om{x). 

For the following let <p(x,y) be a function from C(Ii x V), Y -periodic in y. We will consider 
a(x,y)(j)(x,y). Important particular case: cj>(x,y) = 1. 

Proposition 4.4. 1) If om(x) is measurable in then a(x, x/e)<f>(x, x/e) is measurable in fl. 2/ 
If om (x) G L 1 ^), i/ien a(x, x/e)/(x, x/e) G L 1 (fl). 

Proof. We only have to consider the case e/e 

1) 4>(x, x/e) G C(fl) is measurable. a(x, x/e) is measurable since it is a sum of measurable functions 

©• 
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2) If II is bounded with the diameter 2 R, then U/ej e A/ is bounded with the diameter D := 

2(R + VdA), 


£ d £ d 




/eJ s 


/6J e 


/GJe 


c-d Ip -Id 

E ^ < I LJL D d 


/GJe 


£ d £ 2d 


W(xi) := |z = x/ + -(x — xj) | x G A/j = jz = x/ + -h | x/ + h G A/j C 


i(x,x/e)|| L i (n ) = / |a(x,x/e)| dx = E / |a M (i/ + -(a: - x/))| 

n W, 6 


dx < 


< 


_^ r _^ g-ci 

E^ / - IWIz^n) E? - 11 “a/11^1(0) 


/GJe 


W(xj) 


/GJe 


16 “ e| -Z/. 


g2d 


|a(x, x/ e)<j)(x, xj e)| dx < \\f)\\c / |a(x,x/e)|dx < M\c\\aM\\ L i^) D ° 


Proposition 4.5. If o,m(x) G L 1 (i2) then 

M(x) = f a (x, y)<j>(x, y) dy, M+(x) = ( |a(x, y)4>(x, y)\dy are continuous in Q, 


□ 


IY 


IY 


=d 


M(x),M+(x) are bounded by ||</>||c'||aA/|| L i(n)/' 

Proof. a(x, •) G L^Y") since it was constructed from <zm(■)■ Therefore M{x) and M + (x) are well 
defined. To show continuity let us fix an arbitrary E > 0. 


\M(x + h ) — M(x)| = |y a(x + h, y)(f>(x + h, y) dy - J a(x,y)c/>(x,y)dy 
j a(x + h, y) [<f(x + h, y) - (j>(x, y)\ dy + J a(x + h, y)cf>(x, y) dy - j a(x, y)<j>(x, y) dy 


Y 


Y 


Y 


For continuous 4> one can find such di that |</>(x + /i, y) — </>(x, y)| < i?e d /2||aM||^i(n) when |/i.|co < di 

(|/i|oo = max |hfc| we distinguish from the vector’s absolute value \h\ = \/Ylk h|)- This means that 
fc V 

the first absolute value is less than E/2. 

Now we consider the second absolute value. Using that a(x,y) = a(x,ey), substitution of variables 
z = ey we obtain 

/ a(x,y)(j>(x,y)dy = — } d(x, z)(f>(x, z/e) dz = 

Jy £ JsY 

d(x, z)<f{x, z/e) is e-periodic in z, integral over eY is equal to integral over any e cube 


1 


Similar 


IY 


[ a(x, z)f>(x, z/e) dz = —r f aM(z)(f>(x, z/e) dz 
Jw{x) £ JW{x) 

a(x + h, y)(f>(x, y) dy = —r [ a M (z)(f>(x, z/e) dz 

£ JW(x+h) 
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£ d 


J aM{z)4>(x,z/E)dz — J aM(z)4>{x,z/e)dz 

W (x+h) W(x) 


M 


~ E d 


c 


\a M (z)\ dz 


W{x+h)AW(x) 


fi(W(x + h)AW(x)) < 2d£ d ~ 1 \h\ 00 Using absolute continuity of Lebesgue integral, there exists 
5 2 : \h\oo < <^2 guarantees that the second absolute value is less than E/2 and consequently for 
|/i|oo < minjd'i, 62 } we have \M(x + h) — M(x)\ < E. 


2) \M{x)\ < J \a M (x)\dx = U\\c\\a M \\ L i { ci), 


A 1 . 


Similar with M + (x). □ 

Proposition 4.6. If o,m(x) is measurable in U then a(x,y)<f(x,y) is measurable in Cl x Y. 

Proof. cj>(x, y) is continuous hence measurable. To show measurability of a(x, y) we will construct 
a sequence of measurable functions {ag(x,y)} converging to a(x,y) a.e when 5 —> 0. Let us divide 
into cubes df = i\8, (i\ + 1)5^ x • • • x i^5, (id + 1)5^ , i G T. d . I 5 is a set of indexes i G 7L d 

that df n (1/0. 5 is small enough that He (J df C Cl. For i Gig let xf be an arbitrary point of 

iei s 

df n Q (e.g. the center). 

The function ag(x,y) := a(xf,y), when x G df FI Cl, * G Ig, 

is measurable in fi x Y since a(xf,y) is a measurable function in Y = (0, l) d and df n fi is a 
measurable set. We have to show that the sequence pointwise converges to a(x, y) in U x Y \ O, 
where O := {(x, y) G Cl xY \ y G Af x (x)} is a zero measure set. 

On (df x y) c Of := df X (y n U o c U of. 


of is a measurable set, /xxxv(Of) = /xx(df) x p.y(Y n (J Af x (x)) < nx(0 f)d§ • 

xeinf 


6 S 

Pxxy(O) < ^2 d~nx( Of) < d-na (Cl) -> 0, when 5^0 => Pxxy(O) = 0. 

iei« 

Let (x,y) G (ClxY)\0. It means that dist(y,Af x (x )) > 0 (here dist(y,y) = \y — y\oo)- If we consider 
a 5-partition with 5 < e dist(y,Af x (x)), x G df for some i then (x,y) ^ Of since for all x G df, 
y is enough far from Af x (x). As we know from Prop. POl a(-. v) is piecewise constant in Cl and it 
changes value at those x that y G Af x (x). The whole set df fl Cl belongs to the cube where a(-,y) 
is constant. As a result: Vx G df n Cl, a(x,y ) = a(xf,y). On the other hand from the definition of 
ag : Vx G df n Cl, ag(x,y ) = a(xf,y). Consequently for our particular point (. x,y ) G (df xY)\0 
and small enough 5 we have ag(x,y ) = a(x,y). d 

Lemma 4.1. Let om(x) G A 1 (0); <p(x,y) G C(Cl x Y), Y-periodic in y; a(x,y ) is the Continuous 
Extension of aM(x). Then 


lim 

£—>0 


a(x, x/e)(j)(x, x/e) dx 


Y 


a(x,y)<f>(x,y)dxdy. 


( 10 ) 
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Proof. 1. a(x,y)cj>(x,y) is measurable (Prop. mill . M + (x) is continuous (Prop. 1331) . From Fubini’s 
theorem a(x,y)(f(x,y) £ x Y ), the right hand side is well defined. 

2. Integrals in the left hand side are well defined fProp. WM . 

3. Now we need to check the equality m■ Let E > 0 be an arbitrarily small number, for some e 
(e < e/2) we consider a subdivision of 12 (already defined in Prop. ITdll with ’central’ points xj: 


n= \J Ar, J £ = If* U (j e \ J* nt ); Ji nt = {I e j e I A/ C ny Mr* = y a,. 

Ie Je l£Ji nt 


For not too bad (912 and small enough e, ^(12 \ Q int ) is arbitrarily small: 

~d W 

E A<i < ,|iiiii I, -• 

We will approximate the integrals over hi using the integrals over 12 mt . Let us estimate the errors 
in a similar way as it was done in Prop. m Prop. PI 


ln\n int 


\a(x,x/e)<j)(x,x/e)\dx < \\4>\\c Y rzr / \aM(z)\dz< 

,_rr Jwixj) 


< 


\E — £\ 


E 6 E 


— ll'f’llc||aM|| i i(n) 2d 


5||0||c||aM|| 


i£j s \J int 
< 


LpO) 


E/5, [ \M(x)\dx <E/5. 

Jn\n int 


Since M(x) is continuous in 17, then for small enough e, the right hand side integral in (II HI) can be 
approximated by the sum 

a(x,y)(/)(x,y)dydx^ Y fi(A /) / a(i;j, y)4>(xi, y) dy = 


i& 3% 


Y 


E 

/G Jl™ 


£ d 


[ a(x I ,z)<j>(x I ,^)dz = Y TWpK T [ 


, /m/ * \\ / a M (z)(l)(x I ,^)dz. 

IeJ 1 nt P( w ( x i)) Jw(±j) £ 


( 11 ) 


with error not greater than 2E/5. 

The integral in the left hand side of © can be approximated by 


/ a(x, x/£)(j)(x, x/e) dx « / clm\xi 

Ju. JgJmt A r 


+ - (x - £/ 


x, —), dx = 


or in the new variables 2 = xj + (x — xj)e/e, x = xj + (z — xj)e/e, A/ —> W(xj): 

-d 


^ h l~ 


JeJ mt fc -'W(xi) 

_ m(Aj) f 

j£jint /r(hF(x/)) Jw(xi) 


Xl z - Xl 

x{z), — H--— 


dz = 


aiw(z)</>(x(z), dz. 


( 12 ) 


The last equality is since x/(l/e — 1/e) = I, cj>(x,y) is T-periodic and e/e = /x(A j)//x(W(xj)). The 
approximation error is not greater than i?/5. Further approximation of (II2D : 

M A /) 


E 

/gj™* 


H(W(X!)) Jwixt) 


a M {z)mxi, ^) dz 

r ^' 


(13) 
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has an error 


E 

ieJi nt 




/ a M {z) 

4>(x(z), - cfixr, 

>W(x /) 

\ £J \ £/J 


dz 


which we can estimate in absolute value like in Prop. 14.41 (restricting to e < e/2): 


5||aM 


E p < < E/5, 

jgjint 


where |0(xi,y) - ^(x 2 ,y)| < b = -Ee d /5||aM|| L i(Q)M^) when |xi - x 2 |oo < A/2 < e. 
Now we compare m and m- 


For I G W(xj) = {i 6 I d \x — x/|oo < 

\d =d 


£ 2 


W(xj) C W(X!). 


2|e — e | 

/i(W \ W) = //(VF) — = e 2d /(e — e) d — 0. Then for small enough e 

1 


f) * - -fiw) f w I) dz 


< 


»(W) 

rtw) - n(w) r 
n{W)n{W) Jw 

< (»(W) - y{W)) - Wd -FT 


< 


|dz + 


1 


| a M (z)<t>(xi, dz < 


p(W) Jw\w 

° /L \o>m(z)\ dz < E/5p(Q). 
Jw\w 


^ t{W) Jw\w 

For the second term we used the absolute continuity of Lebesgue integral. 


For small enough s: 

| / a(x, x/e)(j>(x, x/e) dx 

Jn 


n JY 


AE E 

a(x,y)<j)(x,y) dy dx\ < — + - T777 ^ ^ m(Aj) < E. 


5 5/i(12) 


ieJi 


□ 


5 Properties of the two-scale Discrete Extension 

In this section we deal only with the extension a(x, y ) constructed from clm{x ) in the Subsection 12.21 

Proposition 5.1. For fixed x G P n Clj, a(x, ■) was constructed piecewise from aM(')> namely 
is divided into l d cubes by the grid 

Af x (x) = jy G R d | 3k G {1... d}, i G Z : y k = x{/e + i — 1/2 j , 

each cube corresponds to the same e-cube Wj. If x i,x 2 G Qj then Af x {xi) = A/" x (.t 2 ). 

Proposition 5.2. Let us now fix some y G M d . The function a(-,y) is a constant in each ilj. 
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Like in the previous section, we need a simple representation of a(x,x/e). For each Wj if is 
convenient to correspond an e-periodic function 

dj(y) := a M (y) when y G Wj. Periodically expanded to W l . 

a,j(y) := dj(sy) is 1 ^-periodic in M d . 

The simple representation of a(x, y) is 

N N 

a(x,y) = ^2 1 n j (x)a j (y), then a(x,x/e) = ^ 1 n -(x) aj(x/s). (14) 

j =i j =i 

To describe behaviour of aj(x/e) inside Qj let us define the following e-cubes in 

with centers in x\ := x^e/e + Is Aj- = {x G M d \ \x — x\\oa < e/2} 

If x G A} then x = x\ + h (|/j|oo < e/2), 

aj( x/e) = cij(x 3 + el + he/e) = {e-periodicity} = dj(x J + he/e) = o,m(x 3 + he/e) 
since x° + he/e G Wj . We also have a partition of ttj: 


f lj = [J A/, / belongs to J £ (j) when A} := A} n fij / 0. 


We have proved the following 

Proposition 5.3. The simple representation of a(x,x/e) for all e > 0 is : 

N _ 

a(x,x/e) = Y Y 1 A j(x)a M {x j +-(x - x°j)) 
i =1 ieJ e (j) 


We again assume that f(x,y) G x Y), Y -periodic in y. 

Proposition 5.4. l)If om{x) is measurable in then a(x,x/e)c/)(x,x/e) is measurable in Q. 2) If 
om{x) G T 1 (fl) then a(x, x/e)<f)(x, x/e) G L 1 (fl). 


Proof. 1) aj(x/e) are measurable; Qj are measurable sets; a(x,x/e) is a sum of measurable func¬ 
tions. </>{x, x/e) G C(fl) is measurable. 


N r - N 

2 ) / |a(s,s/e)|<^ Y j. . |om(^ + ^{x - xj))| = Y Y 

3= 1 /£Je(i) 




e\ d 


£/ JWi 


\<3m\ < 


< 




14(0) 


AT 


e d 


E E 

3=1 ieJs{j) 


e d = 


\\ a M\\ 


24(0) 


N 


Y^{ U ^i)< N \\ a M\ 


Z4(fi) 


e + 2e 


3 = 1 l€ J £ (i) 

the measures were estimated by (e + 2e) d since fij C Wj and Wj, A} have sides e, e respectively. 
Therefore 


\a(x,x/e)(f(x,x/e)\dx < \\4>\\ c J\a(x,x/e)\dx < A r ||</>||c||aM|| L i(n) 


e + 2e 


□ 
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Proposition 5.5. If o,m{x) E L 1 (0) then 


M(x ) 


a(x,y)4>(x,y)dy, M + (x) 


a{x,y)<t>(x,y)\dy 


JY JY 

are continuous in each f lj and M{x),M + {x) are bounded by ||^||c?||® jw|| x,i(n)/ ^■ 
Proof. 


M(x)= [ a{x,y)<f>(x,y)dy = [ aj(y)(f)(x, y) dy = -j [ aj(z)<l>{x, z/e) dz = 
JY JY £ JsY 


= {dj(z)(J)(x, z/e) is e-periodic in z) 




o-M(z)(j>(x, z/e) dz. 


Let E > 0 be an arbitrarily small number. For X\,X 2 E ttj, \(/)(xi,y)—<j)(x 2 ,y)\ < e> = Ee d /\\aM\\ L i^, 
when \x\ — X 2 I 00 is small enough: 


\M(xi) - M(x 2 )| < 1 / |a M (^)||^(^i,^/£) - </>(a; 2 ,z/e)| dz < E. 

hK w j) JWj 

\M(x)\<^^~ [ \a M (z) \ dz < f \a M {z)\dz. 

£ JWj e Jh 

Similar with M + {x). □ 

Proposition 5.6. If czm(x) is measurable in 0 then a(x,y)(f(x,y) is measurable in Pi x Y. 

Proof. a(x,y) is a sum of measurable functions (fill) . 4>(x,y) is measurable. □ 

Lemma 5.1. Let om(x ) E L 1 (fl),- <f{x,y) E C(Ll x Y), Y-periodic in y; a(x,y) is the Discrete 
Extension of qm^x) ■ Then 


lim / a(x, x/e)<j>{x, x/e) dx = / / a(x,y)<j)(x,y) dx dy. 
£ ^° Jn Jn jy 


(15) 


Proof. Left and right integrals are well defined. Let us fix a small E > 0. 

N 


I / a(x, y)<f>{x , y) dxdy = V" ] / a M (z) / </>{x,z/e)dx dz. (16) 

Q.JY hyWj) JWj l Jn. j 


Let J* nt (j) consists of those / E J e (j) that Aj c Llj. Llj 11 = U /gJmt / N Aj. For not too bad 


boundaries dPlj one can find such small e that 


N N 

E \ sif) = E E 

i =1 i =1 /GJ e (j)\Ji"*(j) 




Ee 




2 W AM 


Li(P) 


Therefore the error of the approximation 

N 


/ a(x, x/e)4>(x, x/e) dx = / , 

isi rTT) ■ 

.7=1 /GJrfl 


’ om{x^ H—(x — Xj))(j)(x, x/e) dx 
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(17) 


N 


j =i ieJi nt (j) ^ 
can be estimated in absolute value as 

N 


^2 ^2 . a M (x j + ~{x - i? I ))<f>(x,x/£)dx 


E E 

j = l/GJeO)\J|™ t (i) 


| c[ \a M (x> + -(x-d? I ))\dx< 

J A] £ 


N 


< 


E E ■ 

0 = 1 ^eJ e (i)\J* nt (i) 


\c d 


J \aM(z)\ 


dz < 


||0||c||aM|| L i ( n) 


N 


£ d 


E E 


Wi 


3 = 1 /GJ e (j')\Jj nt (3) 


e d <* 

~ 2 ' 


We continue EJ by introducing new variables for each integral over Aj: z = V + § (x — xj), 
x = xj(z ) = xj + §(2 — x- 7 ). Additionally we use that x/e = z/e + / and (f> is Y- periodic. (UTI) is 
equal to 

TV 

'£\ d 


N 

E E 

3 = 1 


aM(2)^4(4 7 ) = 


w,- 


N 

E 


- mte 


f a M (z)[ Y ^(^(z),^ 

^ r T -n -f 




2 ; 

£/J 


dz 


and it is approximately equal to dm) . If z G Wj then xj-(z) is some point in e-cube Aj. For small 
enough e the integral from the continuous function f> over iij can be approximated by a sum with 
an error not larger than 5 for all z G Wj . j G 1 ,,N: 


I f <l>(x, -2 )dx- Y £ d <l>(fi(z), |)| < <5 = 

Finally, for small enough e 

/ xx E <5 

a(x, -) 4 >(x, -)dx\ < — + Y 

o 7—1 ' J 




2A||aM|| L i(Q) 


|om(z)| dz < E. 


n y 


Wi 


□ 


6 Admissibility of the Continuous , Viscrete Extensions 

Starting from here, if it is not explicitly mentioned, then a two-scale extension means either the 
Continuous or the Viscrete Extensions as defined in Subsections 12.11 l2~2l 

Corollary 6.1. If a m(’) £ P £ N, f G C(Cl x Y) is Y-periodic, a(x,y ) is either the 

Continuous or the Viscrete Extension of clm{x) then 


lim / a(x, x/£) p cj)(x, x/e) dx 

Jn 


Y 


a(x, y) p (f(x, y) dx dy. 
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Proof. \om\ p = \^m\ then &m(■) = clm(-) p € T^P). As we know from Cor. 12.ll the two-scale 
extension for 6 m( - ) is b(x,y) = a(x,y) p . Then we can apply hem. fCTI or 15 . II for 6 m(-)- □ 

Corollary 6.2. If clm(-) £ T 2 (P), V’ £ C(^ x Y -periodic then a(x,y)'i/j(x,y) is an ’’admis¬ 

sible” test function . 

Proof. We should check (HU, (l8bl) for a(x,y)'i/j(x,y). P is bounded; therefore clm{-) G T^P). For 
arbitrary f> G Btf we can choose p = 1 and (p'f G C(HxF) instead of </> in Cor. 1(111 to verify (Hal) . 
The second condition (I8bl) is again a consequence of Cor. lO.ll with n = 2 and <f = i/> 2 £ C(QxY). □ 


7 Application to the elliptic equation 


We return back to the practical problem from the Subsection 12.31 

In the context of two-scale convergence, the sequence of problems © was investigated in Q],§2. 
Now what is required is to go through the proofs in order to convince ourselves that they still work 
in our case when a(x,y) is a two-scale extension of affix'). 

a|£| 2 < f, T a M (x)f, \a M (x)£\ < /?|f| in P, for any 

implies that (Cor. 12.111 

a|£| 2 < £ T a(x,y)£, \a(x,y)£,\ < /3|£| in P x Y, for any £ G 

For any e > 0, a(x,x/e) is measurable in P provided om( 0 is measurable in P (Prop. Prop. IQ1 
<f> = 1). Therefore for / G L 2 ( P) the problems © are uniquely solvable and their solutions are 
uniformly bounded in Hq(Q). 

Theorem 7.1. The sequence u e of solutions of converges weakly in Hq(Q) (and strongly in 
L 2 ( P)j to no, a unique solution of the limit problem: 


where 


V°: - V • (A(x)Vn 0 ) = f in P, n 0 |an = 0. 

Mx) = f Y V) {yy w j{ x , y) + ej) dy , 


ej - basis vectors, Wj(x,y) (j = 1 ,d) are solutions of the cell problems: 


-V„ 


(a(x,y)(y y Wj(x,y) + ej) j = 0 in Y 

f Wj(x , y) dy = 0, Wj(x, y) is Y-periodic in y 
Y 


(18) 

(19) 


( 20 ) 


Proof. For bounded P we have L°°{ P) C L 2 (P). f>(x) G D{ P), (f>i(x,y) G D[Q-,C™ r (Y)\ have 
continuous derivatives in P x Y and according to Cor. 16.21 

[V0(x) + \7y4> 1 (x,y)} T a(x,y), [V x c(> 1 (x,y)] T a{x,y) 


are row vectors consisting of admissible test functions. Therefore it is still possible in this case to 
pass to the two-scale limit in fl], (2.10) to obtain |Tj, (2.11). The remaining part is given by flj, 
Proof of T.2.3. 

We note that the uniqueness of the solution to the limit problem resulting in the convergence of the 
whole sequence (not just some subsequence) is important to insure that the solution to the initial 
problem V = V c belongs to the convergent sequence. □ 
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What can we say about the averaged coefficient A1 In the case of the Continuous Extension 
a(x, y ) this coefficient should be calculated at each point x G and it depends on the initial 
coefficient aw(-) in e-cube W{x) around x. What happens with A if we slightly move from x to 
x+h? For small enough h the volume W (x) has a large intersection with W (x+h) and consequently 
the coefficients a(x, •), a{x + h, •), which play a crucial role in the cell problem, differ from each 
other only in a small volume. Therefore continuity of the averaged coefficient depends on the form 
of cell problem. 

Proposition 7.1. The coefficient A(x) calculated from the Continuous Extension a(x,y) is con¬ 
tinuous in O. 

Proof. Although our cell problem (needed for calculation of A(x)) is formulated in Y and has a 
variational form: find Wj(x,-) G H^ er {Y) \ R such that 

J^y4>(y) T a(x,y)V y w j (x,y)dy = - J 'V y (/>(y) T a(x,y)e j dy V</> G Hp er (Y) \ R, 

we prefer to deal with the cell problem in terms of om{- ), not in a(-, ■). To do this it is better to 
substitute e-periodic functions for T-periodic: 

~ z z z 

cj)(z) := Wj(x, z) := Wj(x, -), a(x, z) = a(x, —) 

After this substitution the integrals will be over eY from e-periodic functions in z. Therefore they 
are equal to integrals over W(x), where a(x,z) = om{z). 

In new terms the problem has the form: find Wj(x,-) G H} )er (W (x)) \ R such that for all 
4>(z) G Hp er (W(x)) \ R holds the equality B(uij(x, •), <j>) = Cj(4>), where 

B(w,4>):=[ \' z <f)(z) T a M (z)\ / z w(z)dz, Cj{<j>) := [ V z (f{z) T a M (z)ej dz. (21) 

JW(x) £ Jw(x) 

The e-periodic function from H l per (W (x)) \ R is also a function from Hp er (W) \ R, where W is an 
arbitrary e-cube, ||</>||ifi er (^( x ))\R = U\\H^ er (w)\R- Th e Poincare inequality: for all f> G H^ er (Y)\ R 
and f G H^ er (W ) \ R 

/ Hy) 2 dy < C# [ \S7 y (t){y)\ 2 dy, [ f>(z) 2 dz < C|e 2 [ \V Z f)(z)\ 2 dz. 

JY JY JW JW 

We note that the ’small’ size of e is not important here. It is just a fixed constant. The bilinear 
form B is elliptic and bounded on Hp er (W) \ R: 

B{w,w) > a J w ^ \S/ Z w{z)\ 2 dz > l + a c 2 ^ \MH}, er {wmi 

\B(w,f>)\ < P\\w\\ H i er (W)\R\\&\\Hi er (W)\R- 
The linear functional Cj is bounded on Hp er (W) \ R: 

\£j(4>)\ < 1 ||^llHi er (vy)\K- 

Therefore the cell problem has a unique solution wj(x, •) G Hp er (W) \ R, satisfying: 

A + Cle 2 

I \Wj(x, ')ll/fi er (W)\R < --- 
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The formula m written in terms of Wj is 

Aj(x) = X [ efa M (z)(s'V z w j (x,z)+e j )dz. 

£ Jw(x) 


( 22 ) 


x is an arbitrary point from 12. To check the continuity of A(-) we fix some point x £ 12 and consider 
some point x + h € 12 to compare A(x ) and A(x + h). The cell problem for A{x + h ) is 


/ 

Ju 


1 


V z (j){zY clm(z)\ 7 z Wj(x + h, z) dz = —- / V Z (j)(z) 1 a,M{z)ej dz (23) 

W{x+h) e JW(x+h) 

S + := W{x) \ W(x + h ), 5_ := W(x + h)\ W(x). W(x) = W(x + h) U S+ \ 5_. = 

Iw(x+h ) + fs+ ~ Is ' The P r °blem dUD can be re-written: 

I V z (f)(z) T aM(z)'V z Wj(x,z)dz = ~ / 'V z 4>(z) T aM(z)ej dz+ 

JW(x+h ) £ JW(x) 

+ / V z (t>(z) T a M (z)V z Wj(x,z) dz - / V z (j(z) T clm(z)V z Wj(x, z) dz (24) 

./s_ ./s+ 

We substitute (EH) from (EH denoting 0(z) := Wj(x, z) — Wj(x + /r, z): 

[ z ^(z) T a M {z)V z 9{z)dz = (25) 

Jw(x+h ) 

£(</>) = t: [ V z 4>(z) T ciM{z)ej dz[ V z 4>(z) T a,M(z)ej dz+ 

£ JS- e 4S + 

+ / V z ct>(z) T a M {z)V z Wj(x,z)dz- / V z cj)(z) T a M (z)V z Wj(x, z) dz 

Js- Js+ 

The point x was fixed. Consequently the function ibj{x , •) is also a fixed function. 0 belongs to 
i7p e r(lT) \R. Its norm can be estimated by treating 112 5 j) as a variational problem for the unknown 


r [ V z <f>(z) T a M (z)ej dz 
£ Js+ 


< 


m er (w)\m- 


vW 


/ V z (t>(z) T a M (z)V z Wj(x,z)dz \ V Z 4>(z)\\V z Wj(x, z)\dz < 

JS+ JS+ 

< /?||V 2 ^(2:)|| [i 2 (s+)]d ^£ 


| V z Wj(x, z )| 2 cfe 


| V z Wj(x, -)| is a fixed function from L 2 (W{x)). Due to the absolute continuity of the Lebesgue 
integral for arbitrary Eh > 0 one can find such 5 > 0 that the integral from |V 2 u >j(x, -)| 2 over any 
set in W{x) is less than Eh if the set’s measure is less than 6 . The measure of S+, is arbitrarily 
small provided h is small enough. 


L 


V z ^{z) T a M (z)V z Wj(x, z) dz 


< 


m er (w)\m 


V~Eh 
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Similar estimations can be done for the integrals over S-. Therefore ||£|| and consequently 
||01| H 1 er (W)\s. are arbitrarily small for small enough h. 


1 + Cl £ 2 


W d \\H^ r (W)\R < 

This helps us to estimate \Aij(x) — Aij(x + h)\. 


\C\ 


a 


Aij(x ) - Aij(x + h) = ^ 


>W(x+h) 


ej clm{z)eV Z 0(z) dz+ 


+Jd, [ eJa M (z)(£'V z w j (x,z) +ej)dz-f ej a M (z)(eV z Wj(x, z) + ej) dz 
^ J s+ ^ J s~ 

The absolute value of the first term and the integral over S + splitted into two parts 


1 

Cd 


efa M (z)eV Z 0{z) dz\ < /3\\6\\m er (w)\R, 


sd/ 2-1 


e Jw(x+h) 

z# j efa M (z)£V z Wj(x,z)dz\ < ( 3 ^n(S + )E h j £ d ~ 1 , 


'S+ 


£ d 


'S+ 


ej a,M{z)ej dz\ < (3[i(S + )/£ 


=d 


can be made arbitrarily small (due to the small terms ||0||jfi (vy)\R, £?/i, n(S+)) by choosing small 
enough h. Together with similar estimations for S- we have the continuity of A,j(x). □ 


In the case of the Discrete Extension a(x,y), the averaged coefficient A(x) is constant in each 
fife (k = 1,..., N). To determine it one has to solve N cell problems ()21I) . (I22|) with x = x k (centers 
of Wk). This case is realizable in comparison to solving the cell problems at each point in fl. On the 
other hand the averaged coefficient being continuous can be interpolated between a finite number 
of points where it is calculated via cell problems. Here one should be careful since for small e the 
averaged coefficient A{x) is as oscillatory as the initial coefficient Increasing e we expect 

A{x) to become a function with more and more slow variations and in the subdomains of fl where 
the coefficient clm{-) can be classified as ’spatially homogeneous’ it might be close to a constant 
coefficient. 

In Section. El we had a restriction on e from above: e should be small in comparison with the 
typical size of fl. Solving the limit problem numerically with some typical discretization step h 
provides a restriction for e from below: roughly speaking, £ should not be smaller than h. 

Solving numerically the large number of cell problems is a time consuming task, which can be 
done in parallel since cell problems are independent from each other and the limit problem. The 
computational resources can be also saved at least in the following cases: 

• clm{x) has slow variations (for example it can be a constant) in some subdomain Q sv C f2. Then 
inside Q sv there is no need to average. 

• um(x ) is e-periodic in Cl# C Cl and the directions of periodicity coincide with coordinate axes. 
Then the constant averaged coefficient inside Cl# can be calculated by solving only one cell problem. 

Additionally one can also try to combine this with other types of averaging: 

• if the micro coefficient can be classified as statistically homogeneous in some subdomain Cl s h with 
known averaged value A s h or 

• if the averaged coefficient A e d in Cl e d is experimentally determined. 

In these cases one can use the coefficients A s h inside Cl s h and A e d in Cl e d instead of solving cell- 
problems there. 
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Figure 1: a) initial coefficient om(■); b) averaged coefficient v4(-); c) comparison u(x) with uq(x) 

8 Some concluding remarks 

The coefficient a(x,x/e) is often used in homogenization as a generalization of the periodic coef¬ 
ficient a(x/e). In this paper we propose a way to correspond an averaged (limit) problem for the 
initial microscopical problem with non-periodic rapidly oscillated coefficient using the results from 
homogenization together with a special choice of a(x, y ) - the two-scale extension of the initial 
coefficient. The results from homogenization (if not formal) usually require some conditions (like 
its smoothness) on a(x,y). The lack of smoothness e.g. in the two-scale convergence method can 
be partially compensated by the ’’admissibility” of the two-scale extensions, so that e.g. for the 
second order elliptic equation the convergence of u e to uq still holds as in the periodic case. 

To show that this approach can be useful we present here a ID example where clm{x) from 
® and A(x) from (ITFsl) are plotted in Fig. da,b) respectively. To calculate A(-), the Continuous 
extension for e = 0.1 was used. The semi-analytical solutions u(x) (solid line) and uq(x) (dots) 
corresponding to f(x) = —3sin(10x) are compared in Fig. 0c). 

In a 2D test presented in [S] a fine scale reference solution to V is compared with a lif 1 -corrected 
coarse solution to V° (this classical correction is described e.g. in |2{ p.76]). A(x) in V° is calculated 
via the 'Discrete extension of a randomly generated smooth 2D function ««(•)■ 

In consequent publications we are planning to present numerical results for the elliptic problem 
in ID and 2D more systematically together with some other two-scale extensions. 
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